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Abstract 

We investigate the laminar shedding of hairpin vortices in the wake of a truncated square 
cylinder placed in a duct, for Reynolds numbers around the critical threshold of the onset of 
vortex shedding. We single out the formation mechanism of the hairpin vortices by means of 
a detailed analysis of the flow patterns in the steady regime. We show that unlike in previous 
studies of similar structures, the dynamics of the hairpin vortices is entwined with that of the 
counter-rotating pair of streamwise vortices, which we found to be generated in the bottom 
part of the near wake (these are usually referred to as base vortices). In particular, once 
the hairpin structure is released, the base vortices attach to it, forming its legs, so these are 
streamwise, and not spanwise as previously observed in unconfined wakes or behind cylinders 
of lower aspect ratios. We also single out a trail of f2-shaped vortices, generated between 
successive hairpin vortices through a mechanism that is analogous to that active in near-wall 
turbulence. Finally, we show how the dynamics of the structures we identified determine the 
evolution of the drag coefficients and Strouhal numbers when the Reynolds number varies. 

1 Introduction 

Flows past a truncated cylinder are involved in a wide variety of problems dealing with the design 
of e.g. buildings, ship funnels or chimney stacks. The typical values of the Reynolds number Re 
considered in these problems are greater than 10 4 . In contrast, though relevant with probes used in 
i ntrusive measurement metho ds ( Biihler et al . 20081) and heat transfer in electronic circuit boards 
(jNakamura fc Igarashil 120041) . the flow dynamics at lower Re, i.e. up to 10 3 , have been scarcely 



investigated. 

In these regimes, four main structures are generated in the cylinder wake: the horseshoe system, 
the free shear layer stretching from the cylinder free end, those stretching from both cylinder lateral 
faces and a system of trailing vortices. By convenience, we assume that the cylinder is mounted 
on the bottom wall. 

The horseshoe pattern results from the interaction betwee n the bounda r y layers arising r espec- 
tively at the cylinder upstream face and at the bottom wall ( Bakerl ll9T9i : Lin et al. , 2008 ) . It is 



formed by a system of swirls generated by flow separation at the cylinder front that spirals around 
the cylinder. It remains steady fo r Re < 1500 without any restriction on the thickness of the wall 



boundary layer (jLin et all |2008|) 



The system of trailing vortices consists of two pairs of counter-rotating streamwise vortices 
located below the cylinder tip and above the cylinder base, respectively. According to their location, 
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they are referred to as tip and base vortices. In both cases, their origin still remains an open 
question. From experimenta l flow visua lisations in t he wa ke of a truncated circular cylinder at 



Re > 10 4 , lEtzold fe Fiedlerl (|l976h and IPark fc Led (j2000l) suggested that the tip vortices were 



generated above t he upper cylinder face fr om the rolling-up of the lateral ends of the upper free 
shear layer, while iKawamura et al. (1984) interpreted them as resulting from the tilting of the 
lateral free shear layers in the vicinity o f the cylinder f ree en d. From experimental measurements 
in the same configuration at Re > 10 4 , ISumner et al. suggested that the tip vortices had 

nothing to do with the lateral free shear laye rs. Investigating exp erimentally the flow past a 
tr uncated square cy finder for 200 < Re < 10 4 , IWang fe Zhoul (|2009h drew the same conclusions 
aslKawamura et all (1984). On the other han d, li ttle information is av ailable on the base vortices. 

Fiedlerl (|l976( ). ISumner et all (|2004) and lWang fc Zhoul (|2009h a gree that they resu l t from 



Etzold 



the tilting of the lateral free shear layers in the vicinity of the bottom wall. lEtzold fc Fiedlerl (|197(: 
furthermore indicated that the base vortices were initially aligned along the spanwise direction and 
then tilted along the streamwise axis between the mid-span and the free end. The cylinder aspect 
ratio 7, defined as the ratio of the cylinder height h to the characteristic length of its cross- 
section d, dete r mines which of the tip or the base vortices prevail in the wake (pumner et "oil 120041: 
Wang fc Zhoul l2009h . Also, i ncreasing the thic kness of the boundary layer at the bottom wall 



strengthens the base vortices (jWang et all . 120061 ) . 



Since the horseshoe structure remains steady for Re < 1500, any vortex shedding occurring in 
this regime is fed by the free shear layers only and the cylinder aspect ratio 7 determines whether 
the vortex shedding is asymmetric, symmetric or absent. For high values of 7, the vortex shedding 
is govern ed only by the lateral free shear layers and an asymmetric Ka rman-like vortex street is 
observed ( Liu et al . 2005 ; Sakamoto fc Ariel . 1983 ; Wang fc Zhoul . 120091) . For intermediate values, 
the transition to unstea diness leads to a symmetric vortex street formed by hairpin vortices aligned 
on the wake centreline ( Sakamoto fc Ariel 19831; Sakamoto et all. 19871). Finally, small val ues of 7 



prevent any kind of vortex shedding ( Kawamura et al. , 1984 ; Nakamura fc Igarashil I2004J ) . 



The critical values of 7 separating these regimes depend on the shape of the cylinder cross- 
section. The latter affects the shape of the free shear layers as well as the p osition of the wake 
struc t ures with respec t to on e another and thereby the mode of vortex shedding ( Sakamoto fc Ariel 



1983t ISakamoto et all . 119871) . Also, the thickness of the boundary layer at the bottom w a ll has 
a noticeable influence on the determination of the critical values of 7 ( Sakamoto fc Ariel . 



1983; 



Kawamura et all . |1984) 



The formation mechanism of the Karman vortices relies on the alternate roll-up and shedding 
of the lateral free shear layers as in non-truncated cylinder wakes. In contrast, there has not been 
yet any clear agreement on the scenario lead i ng to the symmetric hairpin vortex street. From 
experimental investigations. ISakamoto fc Arid ( 19831 ) suggested that both lateral free shear layers 
joined the top one to form a single entity in the near-wake and as the latter beca me unstable, an 
arch-t ype vortex was formed and released in the wake. This view is supported by Hwang fc Yangl 
( 2004J ) who performed numerical simulations of the flow past a truncated square cyl inder with 7 = 
0.5 for Re < 2000 and by experimental flow visualisations by Wang fc Zhoul ( 20091 ) at Re = 221. 
The la tter authors also inclu ded the tip and base vortices within the arch- vortices released in the 
wake. Yanaoka et ~d\ (|2007l) simulated the flow past a truncated square cylinder featuring 7 = 1 
at Re = 500 and claimed that hairpin vortices were originally vortices that detached from the free 
shear layer stretching from the cylinder free end and then grew into hairpin vortices. Until now 
though, the generation mechanism of hairpin vortices in the wake of a truncated cylinder has never 
been the object of any dedicated study so it remains rather unclear. 

In this article, we investigate the wake of a truncated cylinder of square cross-section using three- 
dimensional direct numerical simulations. Our goal is to clarify the dynamics of the formation of 
hairpin vortices. Our approach relics on a detailed analysis of the flow structures in the steady 
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regime. We restrict our investigations to cylinder wakes confined in a duct of rectangular cross- 
section, in which the inlet velocity profile is a well defined paraboloid flow. Under these conditions, 
the scenario we shall bring forward is new and differs from those inferred in previous works in several 
points. In addition, we have detected a set of ^-shaped vortices forming between two successive 
hairpin vortices. We shall also give the details of the mechanism of their formation. The present 
flow configuration and the numerical set-up are introduced in section [2] In section |3l we identify 
the steady flow patterns and their dynamics throughout the steady flow regime. Section |4] is 
dedicated to the description of the formation and release of hairpin vortices at the onset of vortex 
shedding. The formation of the f2-shaped vortices and the effects of increasing Re within the 
unsteady regime are also detailed in this section. Finally, in section we show how the evolution 
of the flow coefficients is affected by the flow dynamics and, in particular, how the appearance 
of secondary recirculation regions on the cylinder lateral faces noticeably affects the viscous drag 
coefficient. 



2 Configuration and numerical set-up 
2.1 Configuration 

We consider the flow of an incompressible fluid (density p, kinematic viscosity v) past a truncated 
square cylinder in a rectilinear duct of rectangular cross-section as shown in figure [T] A square 
cylinder (width d and height h) is mounted on the duct bottom wall at equal distance from both 
duct side walls. The upstream cylinder face is normal to the streamwise direction taken along the 
x — axis. The cylinder axis is parallel to the z— axis. The origin of the frame of reference is located 
at the centre of the upper cylinder face. The duct height (resp. width) along the z— axis (resp. 
y— axis) is 2a (resp. 2b). The present configuration features h = Ad, 2a = 8d and 2b = lOd, which 
yields a transverse (resp. spanwise) blockage ratio (3 t — d/2b — 0.1 (resp. j3 s = h/2a = 1/2) and a 
cylinder aspect ra tio 7 = h/d = 4. (3 t is low so that it does not have any noticeable effect on the 
horseshoe system ( Ballio et ai . 19981 ). while the effects of the spanwise flow confinement shall be 



discussed throughout this study. 

2.2 Flow equations and numerical set-up 

We study the flow using three-dimensional direct numerical sim ulations with a finite- volume code 



based on the OpcnFOAM 1.4.1 framework (jWeller et all 119981 ). The flow motion is governed by 



the Navier-Stokes equations ([TJ and @, written in their non-dimensional form as: 

Vu = (1) 
^ + (u.V)u = ^V 2 u-V P , (2) 

where u and p are the velocity and the pressure fields, respectively, and Re = Uod/v (Uq is the 
maximum velocity at the inlet). 

The numerical domain is 15d long upstream and 30d long downstream the cylinder. A no- 
slip impermeable boundary condition is imposed at all walls through a homogeneous Dirichlct 
condition and a homogeneous Neumann condition is applied for the velocity at the outlet. This 
geometry ensures that the outlet boundary has no noticeable feedback o nto the cylinder wake as 
it is located more than 25d away from the cylinder ( Sohankar et aZ.I . ll998T >. Preliminary numerical 



computations at a prescribed Re are performed in the duct with the cylinder absent. The resulting 
velocity profile of the fully-established flow in the cross-section is used as the inlet boundary 
condition for the velocity in the three-dimensional simulations with the cylinder present. For the 



V. Dousset & A. Potherat 



Hairpin vortices in the wake of a truncated square cylinder 









J 




u 




■ 




^ X 




— ► 


d 





u 



(b) 



1 



2d 



Figure 1: Configuration of the study: (a) top view; (b) side view. 



pressure field, a homogeneous Neumann condition is imposed at all boundaries but at the outlet 
where a homogeneous Dirichlet condition is applied. The flow equat i ons fH) and ([2]) are solved in a 
segregated way and the PISO algorithm as detailed in lWeller et al\ (|1998l ) is used to deal with the 
pressure-velocity coupling. Second-order discretisation using central difference schemes in space 
and a quadratic backward scheme in time has been used on a non-uniform Cartesian mesh. 

The numerical code has been suc cessfully tested on the two-dimensional flow past a square 
cylinder at Re = 1 50 ijDoolanl . l2009h and by simulatin g the three-d i mens ional flow past a non- 
truncated square cylinder at Re = 200 as described in Muck et al. ( 2000l ). The latter test has 
be en performe d on a non- uniform mesh featuring 1.4 x 10 6 points, i.e. slightly less points than 



Muck et al\ (|2000l) . The non-dimensional minim um grid size was s et to 0.03 so that a better 
resolution of the boundary layer was achieved than in lMiick et "aZl (l2000h . whose computations were 
performed on a fully unifor m mesh. We have o bserved mode-A flow patterns together with dislo- 

~" (<2000h . Our computed values of the Strouhal coefficient 
0.164 and the pressure part of the drag coefficient Cd p = 2Frj/(pU^hd) = 1.70 (/ 



catcd Karman vortices as in iMiick et 
fd/Uo 



St 



is the vortex shedding frequency and Fu the force ob tained by the integ ration of the pressure field 



over the cylinder surface) have recovered the values of lMiick et all (|2000f ) within 1% and 3% respec 



tively, which ensures that the present numerical code addresses the dynamics of three-dimensional 
cylinder wakes within a fully satisfactory accuracy. 

We have then implemented a test to estimate the required number of mesh nodes for the 
simulations with the truncated cylinder wake. The main characteristics of the tested meshes are 
provided in table [T] In the plane normal to the cylinder axis, we have taken over the same mesh 
as that used in the test on the non-truncated cylinder wake which have yielded very satisfactory 
results. We have then varied the number of points along the z— axis and the distance 8 Z between 
the cylinder top face and the grid point nearest to it. We have simulated the flow at Re = 200 over 
1000 turnover times t u — d/Uo and derived the total drag coefficient Co = 2F x /(pU^hd) {F x is 
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Meshes 


Ml 


M2 


M3 


Number of nodes along e z 


32 


64 


96 


s z 


0.06 


0.03 


0.02 


Total number of nodes 


7.4 x 10 5 


1.5 x 10 6 


2.2 x 10 6 


e st = \1 - St(Mi)/St(M3)\ 


0.10 


0.03 


/ 


e cd = \l-C D (Mi)/C D (M3)\ 


3.51 x 10" 3 


3.05 x 10" 3 


/ 



Table 1: Main characteristics of the different meshes and errors in drag coefficient Co and Strouhal 
number St relative to M3 mesh at Re = 200. 6 Z is the non-dimensional distance between the upper 
cylinder face and the grid point nearest to the latter face. All meshes feature 135 and 120 non- 
uniformally distributed nodes along the x- and y-axes, respectively. 



the streamwise component of the total force exerted by the flow on the cylinder) and the Strouhal 
number St in each case. We have found that both the errors in Co and St relative to mesh M3 
decreased with the number of nodes, which shows good convergence. In order to save CPU time 
and keep a reasonable accuracy in our computations, we shall perform all our simulations with M2 
mesh. 

A final validation test has been performed on the configuration of lLiu et~al. ( 2005 ) . The authors 
investigated the flow past a truncated circular cylinder with aspect ratio 7 = 10 between two infinite 
parallel impermeable walls placed at the cylinder bottom end and at 5 cylinder diameters above the 
cylinder free end. We have simulated the flow at Re = 100 using the same geometry, computational 
domain and boundary conditions. For this test, we have designed a structured mesh consisting of 
a polar mesh embedded in a square of two cylinder diameter width and a rectangular Cartesian 
grid covering the rest of the computational domain. The mesh contains about 1.8 x 10 6 nodes with 
a non-dimensional minimum grid size of 0.02 between the cylinder upper face and the grid point 
nearest to it. We have computed the Strouhal number at Re = 100 and found St = 0.1431 which 
represents an error of 1.34% compared to the value of St — 0.145 in lLiu et all (|2005). In addition, 



our simulations have recovered a uni form spanwise di stribution of the Strouhal number along the 
cylinder span in agreement also with iLiu et al\ (|2005l) . 

In summary, our code has been successfully validated on two different configurations of flow past 
a truncated cylinder involving either a transverse or a spanwise flow confinement. In both cases 
the unsteady flow dynamics have been accurately captured by the computation of the Strouhal 
number. We shall now consider the truncated cylinder wake in a rectangular duct in which the 
flow is confined along both transverse and spanwise directions. 



2.3 Vortex identification 

To explain the formation mechanism of hairpin vortices, we need to identify flow structures and, 
in particular, vortices. Establishi ng an objective definitio n of a vortex has been the subject of 
several studies in the literature. Jeong fc Hussainl ( 1995 ) have assessed different approaches of 
identification and tracking of vortical structures. One of them is based on the analysis of the 
eigenvalues of the symmetric tensor S 2 + ft 2 , where S and ft are the respective symmetric and 
antisymmetric part of the velocity gradient tensor Vu. In this approach, a vortex core corresponds 
to a pressure minimum not induced by viscous effects nor unsteady straining. It is defined as a 
connected region with two negative eigenvalues of S 2 + fi 2 . A vortex is therefore detected at a given 
location in the fluid domain if the median eigenvalue, d enoted A2, is locally nega tive. This approach 
is particularly efficient at spotting ring- type vortices (jJeong fc Hussainl . 119951) . Unfortunately, it 
delivers no information on the rotation directions of the vortical structures. For this reason, flow 
patterns will also be characterised by iso-surfaces of vorticity, especially in the unsteady flow 
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Figure 2: Steady flow patterns at Re = 100. Iso-surfaces of A2 = —0.05 (a) and kinematic 
streamlines with horseshoe pattern in pink, trailing vortices in green, streamlines rejoining the 
stagnation points on the upstream cylinder face in cyan, head vortex in red, spanwise vortices in 
blue, black streamlines are only deflected by the cylinder wake without noticeably influencing it: 
three-dimensional (b), side (c) and top (d) views. 



regime. In the stead y one, we will analyse the streamlines, since they match the flow trajectories 



(IChong et qllll99Ch 



3 Steady flow regime 



We have performed nine successive simulations at increasing Re for 10 < Re < 400. Unsteady 
flows have been computed over a total simulation time close to 1000t„. Our computations have 
yielded a steady flow regime for Re < 150 and an unsteady regime for Re > 200. In no case d oes 
the resulting flow feature any closed streamlines, in line with the findings of lHunt et al\ (|1978l ). 



3.1 Flow patterns at Re = 100 

Figure EJa) shows the steady flow patterns at Re = 100 depicted by iso-surfaces of A2 = —0.05, 
i.e. less than 0.2% of the absolute minimum in A2, and by kinematic streamlines in figures [UJb-d). 
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In figures HJb-c) , one sees that the black streamlines that are simply deflected around the wake 
structure separate from those that enter it at saddle point Si located at x = — 1/2 — x u in the 
centre-plane (y = 0). All streamlines other than the black ones enter the wake structure. 



3.1.1 Streamlines originating in the upstream centre-plane 

Pink and green streamlines originate in the upstream centre-plane (y = 0), follow a similar path 
around t he cylinder a nd separate at saddle p oint S2 located at x = 1/2 + Xd in the centre-plane 
(y = 0) ( Hunt et~al\ . 1978t Sau et al. . 2003). Pink lines depict the horseshoe pattern as shown 



in experiments by e.g. ILin et al. I (l2008l)~ The front horseshoe swirl, denoted HS in figure Ufa), 
is observed upstream of the cylinder base at x = — 1. The green set of streamlines impacts the 
bottom duct wall downstream of a nodal p oint o f attachment N a just in front of the upstream 



bottom duct wall downstrea m ot a nodal point 1 attachment l\ a just m iront ot the upstream 
cylinder face, as observed in lHwang fc Yana ( 2004 ). It then recirculates behind the cylinder within 
the horseshoe pattern where it reaches two foci, denoted Fi and F2 in figure EJd), and subsequently 
spirals upwards and eventually tilts along the x— axis to create a pa ir of counter-rotat i ng str eam- 
wise vortices. T his corresponds to the base vortices identified by lEtzold fe Fiedlerl d 19761) and 
Sau et all (120031). T he foc i were also detected in numerical simulations by iHwang fc Yang (|2004r ) 
and lYanaoka et all (|2007l ). but were not linked to the formation of base vortices. At Re = 100, 
these vortices are very weak so that they only induce a weak pressure minimum that is difficult to 
capture with the iso-surfaces of A2. Far upstream, in the centre-plane (y = 0), the border between 
the pink and green sets of streamlines is located at z = hi — 7 and the upper border of the green 
set at z = hi — 7. 

Interestingly, we did not isolate any streamwise tip vortices, whether in the iso-surfaces of A2 
or in kinematic streamlines. There are two probable reasons for this. Firstly, the spanwise flow 
confineme nt imposes a strong spanwise velocity gradient similar to that induced in thick boundary 



layers and IWang et all (|2006f ) showed that the thickening of the bottom boundary layer enhanced 



the base vortices at the expense of the tip ones. Secondly, the present cylinder aspect rat io is rather 



low and no tip vortices hav e been observed in simulations run with a low aspect ratio (|Sau et al 



20031 : lYanaoka et all . 120071 ) 



A third set of streamlines (in cyan) originating within the centre-plane (y = 0) just above the 
green set of streamlines at (/12 — 7 < z < /13 — 7) ends at a line of front stagnation points on the 
upstream cylinder face along the wake centreline. A half-saddle point S3, indicated in figure [21b), 
is detected along this line of stagnation points slightly below the cylinder tip. 



3.1.2 Streamlines rejoining the centre-plane downstream 

Blue and red streamlines rejoin the centre-plane (y = 0) downstream of the cylinder. Blue stream- 
lines originate upstream, below those generating the horseshoe pattern [see figure [2Jc)] but slightly 
off the wake centreline [sec figure (2fd)]. On the one hand, the streamlines located just upstream 
of S2 head upwards and back upstream along the wake centre-plane behind the cylinder and sepa- 
rate into two substreams at half-saddle point 54 on the downstream cylinder face [sec figure [2jc)]. 
The lower substream is deflected by this downstream face outwards the wake along the y-axis and 
eventually carried away by the free stream to form a pair of counter-rotating spanwise vortices. 
The upper substream heads towards the cylinder tip where the free stream takes it away along 
the wake centreline. On the other hand, the blue streamlines located just downstream of S2 head 
upwards and downstream until they reach the free stream where they tilt along the a;— axis. 

A gap between the blue streamlines moving upstream and downstream is observed at the rear 
of the cylinder slightly below the cylinder tip where a transverse vortex is seen [red lines in figure 
IHc-d)]. This head vortex consists of a pair of symmetric transverse vortices located at short 
distance off the wake centreline and connected to a single transverse vortex at the wake centreline. 
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Only the right part of the head vortex, denoted HV, is visible on the A2-iso-surfaces from figure 
[Ha). The streamlines enter this three- vortex structure at the periphery of the symmetric vortices, 
spiral towards the axis and exit at the periphery of the centre vortex onto the wake centreline. 
N ote that the head v ortex is not generated by strea mlines curling f rom t he cylinder free end as 



Hunt et all (|l978l ). ISlaouti fe Gerrardl (|l98lh and IWang fe Zhoul (|200 9|) but from streamline s 



circulating upwards from the cylinder bottom half as in lSau et al\ ([20031 ) and lYanaoka et al\ ([2007[ h 
Throughout this description of the flow topology at the bottom duct wall and cylinder faces at 
Re = 100, we have isolated three nodal points (N a , F\, F2), two saddle points (Si, S2) and two 
half-saddle ones ( 5*3, S4). This topology thus satisfies the mathematical criterion established by 



Hunt et all (|1978[ ) 



Finally, two free shear layers arise from the lateral downstream edges of the cylinder and stretch 
on both sides of the wake (see z— vorticity contours in figure [3]). Due to the presence of the bottom 
duct wall and the free shear layer stretching from the top trailing edge of the cylinder, their 
strcamwisc length is reduced at both their ends and reaches its maximum in the vicinity of the 
base vortices. We insist that all the steady flow patterns detected in the present computations are 
generated exclusively by streamlines engulfed under these lateral free shear layers. 



3.2 Effect of increasing Re within the steady flow regime 

As Re is increased within the steady flow regime, the sets of streamlines rejoining the line of 
stagnation points at the cylinder upper face extends downwards, while the upstream spanwise 
extension of the set of streamlines generating the base vortices (green lines) shrinks and that 
feeding the horseshoe pattern (pink lines) broadens. The increase of Re therefore implies an 
increase of hi and a shrinking of both hi and /13. 

The overall shape of the horseshoe pattern however hardly changes: the main swirl HS [see 
figure [2fa ) j drift s slightly upstream as observed also in iBakeif ( 1979h Hwang fc Yand ( 2004 ) and 
Sau et al\ ([20031 ). but £2 barely moves, i.e. x u increases and xa changes little. The base vortices 
gain in strength and entangle with the spanwise vortices in the vicinity of the foci F\ and F2 , which 
get closer to the rear cylinder face as Re is increased. 

At Re = 150, secondary recirculation regions appear on both lateral and top faces of the 
cylinder. The lateral recirculations are made of stretched and t wisted blue streamline s and the 

The top 



( 20071) 



top one results from flow separation at the cylinder top face as in lYanaoka et 
recirculation deflects the outer stream further up, which, combined with the increased strength of 
the head vortex, results in the latter shifting into a position higher than the cylinder tip in the late 
stages of the steady flow regime. In this regime, the head vortex also moves farther downstream. 
Finally, at Re = 150, the pair of symmetric vortices forming the head vortex merges into the 
centre-plane (y = 0) so that the streamlines are taken away downstream once they come out of 
the eye of this single remaining structure. 



4 Unsteady flow regime 



4.1 Description and formation mechanism of the symmetric vortex street 

at Re = 200 

At Re = 200, unsteadiness appears through a periodic laminar symmetric shedding of hairpin 
vortices aligned o n a single row along the w ake centreline. This vortex shedding has been observed 
in experiments by Sakamo to fc ArieTl 19831) a nd Sakamoto et al. (1987) and numerical simulations 
by I Hwang fc Yang ( 2004 ) and Yanaoka et al. ( 20071) . The formation and release of hairpin vortices 
depicted by iso-surfaces of vorticity field and of A2 = —0.5 (0.5% of its absolute minimum) are 
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Figure 3: Vortex street at Re = 200 depicted by iso-surfaces of non-dimensional vorticity ui (a-c) 
and A 2 = —0.5 (d-f) at t = 1349i„ (top), 1351t u (middle) and 1353i u (bottom). uj x = 2.5 (resp. 
u x = —2.5) in yellow (resp. cyan); uj y = 5 in green; lo z = 5 (resp. lu z = —5) in red (resp. blue). 
The red lateral free shear layer is not shown in (b). TV, LV and RV are vortical structures shed 
from the top, left and right free shear layers, respectively. BR and BL are the base vortices. SR and 
SL are the secondary streamwise vortices clearly visible on A2-iso-surfaces but not on vorticity iso- 
surfaces where their locations are indicated by dashed arrows. Red and blue arrows on (d) indicate 
the jets induced by the base and hairpin vortices respectively, whose interaction participates to the 
generation of il-shapcd vortices. Supplementary movie 1 (available at journals.cambridge.org/flm) 
shows an animation of the vortex street at Re = 200. 
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represented at three successive time instants during the formation and release of the hairpin vortices 
at Re = 200 in figures EI a-f). 

At t — 1349t u , one observes a hairpin vortex in the early stage of its formation in the near 
cylinder wake and a freshly released, fully formed, hairpin vortex in the far wake. Three distinct 
vortical structures denoted TV, LV and RV, respectively originating from the top, left and right 
free shear layers, are singled out in figure EJd). TV indirectly results from the mutual interaction 
of the base vortices BR and BL, shown in figure EJa), that generates an upwards jet between 
them and lifts the tail of the top free shear l ayer. This results in an inversion of the curvature 
of the latter which becomes unstably curved (jLioul . 11994 ). and leads to the breakaway of its tail 
(TV). At t = 135 li u , two secondary counter-rotating streamwise vortices, denoted SR and SL 
in figures EJb-f), then rise just upstream of the structure shed from the top free shear layer and 
wrap around the base vortices, while TV, LV and RV have gathered into a single bow of vorticity. 
The subsequent pairing between both the base and the secondary streamwise vortices eventually 
triggers the symmetric shedding of structures LV and RV from both lateral free shear layers [blue 
and red iso-surfaces in figureEIc)] observed at t = 1351i u . At this moment, the head of the hairpin 
vortex is completely formed and taken away by the free stream, while the base vortices detach to 
join the head of the hairpin vortex and form its legs visible on the hairpin vortex present in the 
far wake in figure [3jd) . 

Our computations show that the hairpin v ortices are not genera ted o nly by the destabiliza - 
tion of the top free shear layer as suggested in Hwang fc Yang ( 2004 ) and lYanaoka et all ( 2007 ) . 
Instead, they result from the smooth assembling of structures shed from the initially steady flow 
patterns. The head of the hairpin vortex is formed from the smooth assembling of vortical struc- 
tures detached from the top and lateral free shear layers. The hairpin vortic es observed in the 
presen t sim ulations differ f r om th e arch- type vortices detected in experiments by ISakamoto fc Arie 
(jl983h and IWang fc Zhoul (l2009h . firstly, due to the spanwise velocity gradient imposed by our 
flow confinement, the hairpin head is located at the downstream end of the structure. The free 
stream velocity is higher in the vicinity of the cylinder tip so that the hairpin head, located in 
this region, is carried away further downstream than the rest of the hairpin structure. Secondly, 
the hairpin legs arc almost parallel to the streamwise axis and not aligned along the spanwise 
direction as those of the arch-type vortices. In our configuration, the hairpin legs are indeed the 
base vortices which are oriented along the streamwise axis, whereas those of the arch-type vortices 
are the spanwise vortices shed from the lateral free shear layers. Thirdly, in our configuration, 
neither the legs nor the spanwise vortical part of the hair pin vortex grow a pos teriori from the 
structure shed from the top free shear layer, as claimed in lYanaoka et al\ ( 2007 ). Instead, both 
these structures are present in the very first steps of the formation of the hairpin vortices. This 
shows that the base vor tices play an active part in the vortex sh edding mechan ism as in the flow 



past a tra pezoidal tab (iDong fc Meng , 2004 ; Yang et all 1 2 lh , but unlike in ISakamoto fc Arie 
(|l983h and I Wang fc Zhoul (120091 )7 1 n the last two studies though, the inlet velocity is almost span- 
wise invariant and the shedding results from the interaction between the top and lateral free shear 
layers only. 



4.2 Formation and release of secondary Q-shaped vortices 

We have described how the hairpin vortices arc generated and released in the cylinder wake. In 
the unsteady regime, the structures that take part in their formation were detected by plotting 
iso-surfaces of vorticity and A2 for reasonably high values of these quantities. We shall now refine 
our analysis and track weaker structures by bringing the threshold values of A2 and of the vorticity 
closer to zero for the flow at Re = 200. 

Figure Ufa) shows the three-dimensional vortex street at t = 1353i„ depicted by iso-surfaces 
of A2 = —0.03 and two-dimensional vorticity iso-surfaces in three different planes are presented in 
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Figure 4: (a) Iso-surfaces of A2 = —0.03 at Re = 200 and t = 1353. Fields of non-dimensional (b) 
x-vorticity in plane A (x = 12.9) with \tu x \ < 2, (c) z— vorticity in plane B (z = 0) with \oj z \ < 1.5 
and (d) y— vorticity in plane C (y = 0) with \tu y \ < 1.8 (minimum and maximum values in blue and 
red, respectively). Thick white lines on (b-d) are the intersected A2— iso-surfaces. Hi and H2 are 
two successive hairpin vortices, fii, ^2 and ^3 three f2-shaped vortices, SR and SL the respective 
left and right secondary streamwise vortices and BR and BL the respective left and right base 
vortices (see also figure [3]). 
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figures [DJb-d) . One observes a chain of Sl-shapcd vortices between two successive hairpin vortices, 
denoted Hi and H^. Three such f2-shaped vortices, denoted fii, f^2 and ^3, are detected in figure 
Eta). The strength of these vortices decreases with their distance to the downstream hairpin vortex. 
They are located in (y, z)-planes and their bottom ends are connected to the secondary streamwise 
vortices, denoted SR and SL in figure Eld). The fi-shaped vortices are mirror images of the hairpin 
vortex as they all rotate in the sense opposite to it, as shown by figures HJb-d) . Vortex streets 
involving hairpin vortices and secondary reverse vortices (i.e. rotating in the opposite sense to 
that of the hairpi n head) have been observed in the flow past a tr apezoidal tab in experiments 
( Yang et all . [?()() 1 : and numerical simulations ( Dong fc Mend . 2004 ) . Reverse vortices detected in 



the latter configuration were however located below the line of hairpin vortices by co ntrast with our 
observ ations where hairpin and reverse vortices are all aligned. The simulations bv lDong fc Mens 



(|2004j ) have yet clearly established the il-shape of the reverse vortices as in our case. 

The formation of Sl-shaped vortices follow s from a me c hanism encountered in turbulent bound- 
ary layers in channel flows described in e.g. Zhou et al. ( 19991) . During its downstream motion, 



the head of the hairpin vortex broadens and induces a backwards streamwise jet through the head 
[blue arrow in figure [3jd)]. The pair of streamwise vortices forming the legs of the hairpin vortex 
generates an upwards jet between them [red arrow in figure EJd)]. The shearing between the latter 
jet and the one induced by the head of the hairpin vortex causes the formation of a bridging shear 
layer between both streamwise vortices. This shear layer subsequently rolls up into an f2-shaped 
vortex in the (y, z)-plane. This mechanism is active all along the legs of the hairpin vortex and 
leads to the formation of several fi-shaped vortices between two successive hairpin vortices. 

Importantly, the generation of f2-shapcd vortices is a direct consequence of the streamwise 
orientation of the hairpin le gs. This mechanism is thus absent in configur ations where arch-type 
vortices have been observed ([Sakamoto fc Ariel Il983t IWang fc Zhoull2009h . 



4.3 Effect of increasing Re on the vortex street 

At Re = 250, the mechanism of formation and release of hairpin vortices described previously 
remains unchanged, but the shape of the hairpin vortices becomes irregular and weak spanwise 
oscillations of the base vortices are observed. The latter oscillations disturb the shedding of the 
top and lateral free shear layers so that the vortex street turns slightly asymmetric as shown in 
figure [Ua). For Re > 300, these oscillations gain in intensity. The legs of the hairpin vortices 
do not split at the same z— coordinate anymore, while their head turns into a chaotic aggregate 
of structures shed from the free shear layers. Hairpin vortices are therefore released alternately 
on either side of the wake centreline and so a fully asymmetric vortex street is observed as seen 
in figure [5] (b). The chain of Sl-shaped vortices is still observed in this regime, but their shape is 
irregular too and their orientation follows that of the hairpin vortex that creates them. A sim ilar 
evolution of the vortex street with increasing Re has been described bv lYamada et al. (1993) in 



experiments of the flow past a square plate. The gradual change from a regular to a chaotic vortex 
street for 200 < Re < 400 is illustrated in figure [6] by the respective time histories of the total drag 
coefficient Cd, introduced in section l2~2l 



5 Flow coefficients 
5.1 Drag coefficients 

The variations of the total drag coefficient Cd with Re are presented in figure [T^a). One observes 
that Cd decreases within the steady flow regime and slightly increases in the unsteady one. These 
variations reflect the dynamics of the vortex formation region as in the two-dimensional square 
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Figure 5: Vortex street at Re = 250 and t = U10t u (a) and Re = 300 and t = 1044t M (b). 
(left) Iso-surfaccs of A2 = -0.5 (a) and -0.7 (b) and (right) fields of uj x = ±3.75, ui y = 6.25 and 
uj z = ±6.25 (a) and u x = ±4.25, ui y = 9.25 and uj z = ±9.25 (b) in yellow and cyan, green, red and 
blue, respectively. 
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Figure 6: Time histories of the total drag coefficient Cd at i?e=200, 250, 300 and 400. 
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cylinder wake (jSharma fc Eswaranl 120041) . In steady regimes, it encompasses the pair of spanwise 
recirculation regions, while in unsteady ones it includes the formation region of the hairpin vortices. 
In both cases, this region induces an adverse pressure gradient. The increase of Re within the steady 
regime results in the streamwise elongation of the symmetric steady recirculation regions. This 
enhances the adverse pressure gradient which gradually counter-balances the streamwise pressure 
gradient induced by the duct flow and thus leads to the decrease in the total drag coefficient Co- 
in the unsteady flow regime, the mechanism is reversed as the formation region of the hairpin 
vortices shrinks, hence the increase in Co- Remarkably, the present values of Co obtained with a 
cylinder spanning over half the duct height are about half of those obtained with a non-truncated 



cylinder (jSohankar et all 1999). 



The viscous component of the drag coefficient is defined as Cov = 2F V / '{pU^dh) where F v is 
the viscous force acting on the cylinder. We have reported its variations over the cylinder span at 
different Re in figure [TJb). One can identify three regions alo ng the span: two narrow ones close to 
the respective cylinder ends and a larger one between them (jLiu et al ., 2005). 

Cov exhibits peak 

values just where the flow around the cylinder is fast, i.e. below the cylinder free end and at short 
distance of the bottom duct wall within the horseshoe region. In contrast, Cov remains rather 
constant in the mid-span region where the flow recirculates much more slowly. The corresponding 
values are therefore significantly smaller than the peak values reached in both end regions. One 
observes also that Cov turns negative for (—0.4 < z/h < —0.13) at Re = 200 and that the 
z— range with negative Cov widens as Re is increased. This feature results from the appearance 
and subsequent widening of the secondary recirculation regions at the cylinder lateral faces for 
Re > 150. The associated counter-flow induces advert friction at the side walls that reduces the 
overall value of Cov to th e point of making it negativ e. This effect was detected in two-dimensional 
square cylinder wakes by ISharma fc Eswaranl (]2004f ) . 



5.2 Spanwise lift coefficient 

We define the spanwise lift coefficient C z = 2F Z / (pU^d 2 ) where F z is the z— component of the force 
exerted by the flow on the cylinder. The evolution of C z computed only on the upper cylinder face 
versus Re is reported in figure [T^c) . By definition, F z results from the integration of the pressure 
force on the cylinder upper face and is defined up to a constant determined by the value of the 
reference pressure. As a consequence, only the variations of C z reflect the flow dynamics, while 
its absolute values shall be regarded relatively to the reference pressure. In particular, the sign of 
C z bears no significance. C z decreases in the steady regime and slightly increases in the unsteady 
one. As for the drag coefficient, the decrease of C z results from the lengthening of the spanwise 
recirculation regions at the back of the cylinder. The reason for its increase is the appearance of the 
secondary recirculation over the cylinder upper face that induces an adverse pressure gradient. The 
dcstabilization of the top free shear layer in the unsteady regime for Re > 200 slightly enhances 
this trend throughout the unsteady regime. 



5.3 Strouhal numbers 

The frequency of the symmetric {resp. asymmetric) mode was obtained from the time history of the 
total drag coefficient Co {resp. lift coefficient Cl = 2F y / pU§dh where F y is the y— component of 
the force exerted by the flow on the cylinder). We were then able to calculate the Strouhal numbers 
associated to the respective frequencies of the symmetric and asymmetric modes. Their dependence 
on Re is presented in figure \7^d) . At Re = 200, the vortex street is indeed perfectly symmetric 
so the asymmetric mode is absent. At Re = 250, a slight asymmetry is observed in the wake 
due to the appearance of the small vertical oscillations of the base vortices. These oscillations are 
enhanced when Re is further increased until the wake becomes chaotic and completely asymmetric. 
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Figure 7: (a) Total drag coefficient Co on the cylinder stem versus Re and (b) spanwisc distribution 
of the viscous drag coefficient Cdv at different Re (< Cdv >z is the z— average value and the time- 
average value of Cdv is used in the unsteady flow regime) . (c) Spanwisc lift coefficient C z on the 
cylinder upper face and (d) Strouhal numbers St associated to the symmetric and asymmetric 
modes versus Re. The dashed line on (a) and (c) separates the steady and the unsteady flow 
regimes. 

When the wake is symmetric at Re = 200, only the symmetric mode is present with an associated 
St = 0.07. When the asymmetric mode appears at Re = 250, its associated St is much lower 
(St = 0.01) but increases with Re. As the wake becomes more asymmetric, the Strouhal number 
associated to the symmetric mode collapses down to values below those of the asymmetric one. 

6 Conclusions 

We have investigated the flow past a truncated square cylinder in a duct of rectangular cross- 
section using three-dimensional direct numerical simulations for Re < 400. Figure |S] summarises 
the main flow regimes we have encountered. The steady regime has been thoroughly scrutinised 
and the flow patterns identified. The present computations have shown that these patterns were 
all generated by sets of streamlines flowing around the cylinder foot and underneath the lateral 
free shear layers. In particular, in our configuration, the spanwise flow confinement induces an 
important spanwise shear that promotes the development of the base vortices at the expense of 
the tip ones. Also, we have singled out the appearance of secondary recirculation regions at the 
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Figure 8: Summary of the main flow features observed in the present simulations. 



lateral and upper cylinder faces for Re > 150. 

At Re = 200, the flow unsteadiness leads to a symmetric vortex street formed by a single row 
of hairpin vortices aligned on the wake centreline. We have detailed a formation mechanism of the 
hairpin vortices resulting from the smooth assembling of structures shed from the initially steady 
flow p atterns. This scenario has never been exposed so far and differs from t hose given in previou s 
works (|Hwang k Yand . 120041 ISakamoto k Arid . Il983t lYanaoka etHH l2007t I Wang k Zhoul . l2009h . 
To partly explain this difference, we have stressed the effects of the vertical flow confinement 
which enhances the base vortices and induces a more efficient entrainment of the head of the 
hairpin vortices. The latter is therefore located at the downstream end of the hairpin and the legs 
are almost parallel to the streamwise axis. Also, we have detected the formation of a chain of 
17-shaped vortices as a consequence of the streamwise orientation of legs of the hairpin vortices. 
The fi-shaped vortices are well known in turbulent boundary layers and their formation stems from 
the same mechanism as in this classic case. For Re > 250, the appearance and enhancement of 
vertical oscillations of the base vortices causes the vortex street to turn asymmetric and irregular, 
while hairpin vortices become chaotic aggregates of structures. 

We have reported the evolutions of a set of flow coefficients with Re. Firstly, the spanwise lift 
coefficient has been found to increase when the recirculation appeared on the top face. Secondly, 
the transition to asymmetry in the unsteady wake has been shown to induce a collapse in the 
Strouhal number that reflected the disappearance of the symmetric mode associated with the 
hairpin vortex shedding we have identified. 

Finally, although we have clearly identified the dynamics of the hairpin vortices, questions 
about the respective influences of the various parameters set in our numerical computations are still 
open. Firstly, the influence of the cylinder aspect ratio 7 would require a parametric investigation. 
Although it is well accepted that 7 determines whether the onset of vortex shedding gives rise 
to hairpin or Karman vortices, the effect of 7 within any of these regimes is yet to be clarified. 
Secondly, we have used physical walls at the outer boundaries of our domain. Though physically 
realistic, both the transverse and spanwise confinements should be subject to a thorough scrutiny 
to determine their effects on the stability and the shape of the hairpin vortices. The transverse 
confinement in the present configuration is low enough to neglect its effects. One may however 
expect valuable information on whether moderate transverse blockage improves the stability of 
the hairpin vortices. By contrast, we have shown that the spanwise blockage had a significant 
influence on the flow dynamics. One could expect a lighter spanwise confinement to minimise the 
role of the base vortices in these dynamics and result in a modification of the overall shape of 
the hairpin vortices wh ich c ould be closer to the arch-type vortices observed in experiments by 
Sakamoto k Arid (|l983l) and I Wang k Zhoul (|2009h . 



Supplementary movie available at journals.cambridge.org/flm 
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